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ABSTRACT 

Context. In recent years, there has been a number of detections of gradients in the radial velocity profile across jets from young stars. 
The significance of these results is considerable. They may be interpreted as a signature of jet rotation about its symmetry axis, thereby 
representing the only existing observational indications supporting the theory that jets extract angular momentum from star-disk sys- 
tems. However, the possibility that we are indeed observing jet rotation in pre-main sequence systems is undergoing active debate. 
Aims. To test the validity of a rotation argument, we must extend the survey to a larger sample, including younger sources. 
Methods. We present the latest results of a radial velocity analysis on jets from Class and I sources, using high resolution data from 
the infrared spectrograph GNIRS on GEMINI South. We obtained infrared spectra protostellar jets HH34, HH 111-H, HH212 NK1 
and SKI. 

Results. The [Fe II] emission was unresolved in all cases and so Doppler shifts across the jet width could not be accessed. The H2 
emission was resolved in all cases except HH 34. Doppler profiles across the molecular emission were obtained, and gradients in 
radial velocity of typically 3 km s -1 identified. 

Conclusions. Agreement with previous studies implies they may be interpreted as jet rotation, leading to toroidal velocity and angular 
momentum flux estimates of 1.5 km s~' and lxlO~ 5 M yr -1 AUkm s -1 respectively. However, caution is needed. For example, emis- 
sion is asymmetric across the jets from HH 212 suggesting a more complex interpretation is warranted. Furthermore, observations for 
HH212 and HH 111-H are conducted far from the source implying external influences are more likely to confuse the intrinsic flow 
kinematics. These observations demonstrate the difficulty of conducting this study from the ground, and highlight the necessity for 
high angular resolution via adaptive optics or space-based facilities. 

Key words. ISM: jets and outflows — Stars: formation, pre-main sequence, mass-loss, low-mass — stars: individual: HH 111, HH 
34, HH 212 — Infrared: ISM 



[ 1. Introduction of the outflow from Class source HH212, in three slit posi- 

t— I ■ tions stepped parallel to the flow direction, and revealed a dif- 

t— I 1 One of the fundamental problems of star formation theory is f ere nce in radial velocity of 210ns- 1 across the flow in the H 2 

O finding an explanation of how the excess angular momentum 2.12 ^m near-infrared ( near IR) line at 6" (2500AU) from the 

. may be extracted from the accreting gas and dust so that it star dDavis et alll200flh . This was taken as the first observational 

L; ■ may continue to travel inwards and eventually accumulate on the hint f j et rota tion, given the context of a concurrent study of 

. , £h ! newly forming star. With the realisation that observed shocked the HH 212 disk which revealed a r adial velocity gradient in the 

^ . gas was in fact the result of imp acting high velocity jets and same sense dWiseman et all l200lt) . Independantly, other stud- 

H i outflo ws ejected from such stars (|Bally et al.j |2007t |Ray et al.j j es examined jets from less embedded Class II sources at opti- 

. 12007|), came the proposal that such ejecta could be the under- cal and near-ultraviolet (near UV) wavelengths, and harnessed 

lying vehicle for angular momentum transport. The theoretical the high reso lution of the Hubble Space Telescope (HST) in or- 

foundations have long been laid in which jets a r e ejected via der to observe jets closer to their ejection point. The DGTau 

magn eto-centrifugal forces (|Pudritz et al.j |2007t |Shang et all an d RW Aur jets, examined with a similar parallel slit configu- 

|2003), but without observational verification. It became clear ratiorii reV ealed radi al velocity gradients of 5-20 km s' 1 withi n 

that high angular resolution observations are required to test the 10 AU of the star (Bac ciotti et all 12001 IWoitas et all |2005). 

theory, since all such models describe acceleration and collima- Furthermore, these gradients were sustained for ~ 100 AU along 

tion on scales of tens of AU. both i obes G f the RW Aur bipolar jet. This sustainability does 

A significant observational breakthrough has been made in not favour alternative interpretations such as jet precession or 

recent years by our team. We have conducted a series of stud- asymmetric shocking. Such positive indications of jet rotation 

ies revealing for the first time indications that we can observe required a survey. We therefore examined several TTauri sys- 

the jet rotating. This combined endeavour may prove to be the tems, to study a series of 8 jet targets including two bipolar jets, 

long-awaited observational backing for the magneto-centrifugal This time the slit was placed perpendicular to the flow, which 
theory of star formation. The first study constituted observations 
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maximises the opportunity of detecting a rotation signature by 
observing across the full jet width. It also maximises survey ef- 
ficiency and_avoids difficulties introduced by uneven slit illumi- 
nation (IWoitas et al.Ll2005h . Analysis in the optical and near UV 
consistently produced systematic radial velocities differences of 
typica lly 15-25 km s _1 close to the ejection point (ICoffev et al.L 
2004; 2007). These very encouraging statistics led to a pilot 
study to detect gradients in younger Class I sources, since ro- 
tation signatures should be present at all evolutionary stages. 
Indeed, near IR spectra of HH 26 and HH 72 revealed gradients 
in their radial velocity profile s transverse to the flow direction at 
1 000 AU from the star (IChrvsostomou et al.Ll2008h . 

In response to the positive outcome of this pilot study, we 
now conduct a wider survey of younger sources to consolidate 
the results. Three Class 0/1 systems were observed in near IR 
lines, using the high spectral resolution (R~18 000) available 
with GNIRS on GEMINI South. Based on previous detections 
of outflows in the near IR, we choose HH 34 (IDavis et all [200 1 ; 
iDavis et alll2003l HH 1 1 1-H (iDavis etafll2001bl) . andHH 212 
NK1 and SKI dDavis et aUl2000l) . The fFe III 1.64/im emission 
traces the hot inner parts of the jet while the H2 emission at 
2.12/mi traces the warmer outer regions. Analysis of the spa- 
tially broad emission, which is resolved under good seeing con- 
ditions, takes full advantage of the high spectral resolution set- 
ting. Hence, using spectral analysis techniques, we are able to 
determine whether or not a Doppler gradient is present and in 
what sense, and thus the magnitude of the implied toroidal ve- 
locity. 



2. Sample 



2.1. HH34 



The HH34 jet, first reported by Reipurth (1986), is the name 
given to the most recent ejection of the highly collimated S- 
shaped parsec-scale HH 34 flow fro m a Class I source in Orion, 
which is located at 414pc distant dMenten et all 120071) . It has 
since been well studied via imaging and spectroscopy, but per- 
haps most i mpressive are the de ep high resolution HST opti- 
cal images (Reipur th et all 120021) . used to determine morphol- 
ogy, photometry and excitation of the flow, and reveal that 
the jet abruptly changes direction possibly due to the power- 
ful tidal effects of a companion star. This is in addi tion to the 
large s cale S- shape of the flow a lready reported by Ball v et aTl 
(1994j) and iDevine et al.l d!997l) . Periodic time- variable ejec- 
tion modeling of the jet hav e successful l y repr oduced the flow 



structure and morpho logy (Ra sa et ail [1998; 2002c), while 



Masciadri et al.l (120021) investigate precession of the flow as a 
deceleration mechanism. More recently, hi gh resolution opti- 
cal integral field spectroscopy observations dBeck et all [2007) 
show agreement with the kinematics and electron density struc- 
ture predicted by existing internal working surface models, al- 
though radial velocity studies show no evidence of a Doppler 
gradient across the flow but this may be because the effec- 
tive resolution of 20kms~' did not allow a detection. Close 
to the source, the H? counter-part to the optical jet emission 
was first traced by IDavis et alT d2001l) . wh ile the same regio n 
was also seen in [Fell] 1.64/im emission dDavis et all [2003). 
Further recent near- IR kinematic and d ia gnostic studie s of th e 
jet physics include iPodio et al.l (120061); iTakamietall (|2006); 
iGarcia Lopez et al.l d2008l) :lAntoniucci et al.l (l2008lK 



2.2. HH 111 - H 

The highl y collimated HH 111 Class I outflow was first re- 
ported by iReipurthl d!989l) . Located in Orion, the HH 1 1 1 
flow has been studied in detail over the years. Observations of 
the o utflow itself at high resolution include H ST optical im- 
ages dHartigan et all 1200 U iReipurfh et"afl 1 19971) and spectra 
dRaga et all 120021) . As with HH 34, it shows a chain of well 
aligned knots ending in a bowshock, the H knot being one of 
the brightest and located at 33" (13700AU) from the source. 
Originating from a known multiple system, it is remarkable for 
its stability over parsec s cales. It is also a strong H2 emitter 
dGredel & Reipurthl 1 1993b. a nd is a ssociated with a powerful 
CO outflow (jCernicharo et all 1 19961). The H? flow has b een ex- 
amined using echelle spectroscopy dDavis et all 1200 lbl) which 
revealed knot H to have a double-peaked profile, interpreted 
with a simple, geometrical bow shock m odel. A sister bipolar 
flow, HH 121 dGredel & Reipurthl, 1 19931; 1994), originates from 
the same position but is offset in RA. by 60° from HH111. 
The core of this quadr upolar outflow has also been examined 
dRodriguez et all l2008) with a detection possibly indicating two 
disks. Models of the ejection history of HH 111 western lobe 
dRaga et all 120021 iMasciadri et alll2002bl) show the bow shocks 
to be the result of an ejection ve locity time-v a riabi lity, while near 
IR spectral diagnostic studies dNisini et all 120021) demonstrate 
the nature of the shocked gas. 

2.3. HH212 NK1 & SK1 

HH212 is a well-known H2 outflow from a Class source 
in Orion. Early imaging observations dZinnecker et all Il998l) 
clearly show a highly symmetric bipolar flow, with a total 
length of 240"(0.5pc). The H2 knots closest to the source in 
each lobe (~6") are named NK1 and SKI. The H H212 sys- 



tem i s also associated with a collimated CO outflow (Lee et al 
2006|), and an SiO flow ( ICodella et alll2007tlCabrit et alll2007 



Lee et all 2007t 2008). Ne ar IR diagnostic studies include 



Caratti o Garatt i et al.l (2006), wh ich examines physical proper 



ties and cooling mechanisms, and Smith et al.l (I2007I) . which re- 
veals excitation properties consistent with outward-moving bow 
shocks close to the plane of the sky. 

The first hint that jet rotation is detectable in a y oung stellar 
j et act ually came with the observations of HH 212 bv lDavis et al.l 
(2000), based on a difference in the radial velocity across the re- 
ceding lobe of the H2 flow. This was interpreted as such given the 
agreement with the sense of a radial veloc i ty grad i ent across the 
disk o f the same system dWiseman et all l200ll) . ICod ella et al. 

(2007) report no sign of jet rotation in the SiO emission near 
NK1, whereas in the region SK2-SK4 (-10-14" from source), 
they fin d a gradient in a d irection contrary to that reported for 
SKI bv lDavis etafl (120001) . It is therefore also contrary to that 
of the disk which has a 5kms _I difference between the blue- 
shifted ammonia g as in the north-east an d the red-s hifted gas 
in the south-west dWiseman et all 1200 lb . However, iLee et all 

(2008) recently conducted a higher resolution rotation study of 
the HH212 flow and report matching gradients in the southern 
and northern lobes of the SiO jet in a sense that also matches 
that of the disk. They note that the slope of the gradient in the 
northern lobe (at 800 AU from the disk plane) is smaller than 
the southern lobe (at 450 AU). Indeed, a very comprehensive 
H2 study of this sys tem was conducted using GEMINI/Phoenix 
(ICorreia et all 120091) . in which the various kinematic possibili- 
ties are modeled with the conclusion that the observations can- 
not be reproduced by jet rotation alone, though jet rotation may 
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be included with other effect s such as precess ion. Furthermore, 
near IR diagnostic studies (ISmith et all 120071) report a gradient 
in excitation transverse to the jet axis across in the inner knots, 
suggesting a transverse source motion rather than precession, 
possibly related to the jet bending and the transverse gradient 
in radial velocity. A similar gradient in excitation transverse to 
the jet axis in the T Tauri star Th 28 was reported by Coffe v~et~afl 
(2008), which possibly suggests asymmetric shocking. 



3. Observations 

Observations were conducted, in queue mode, with the GEMINI 
Near InfraRed Spectrograph (GNIRS) on GEMINI South in 
mid-October 2006, late December 2006 and early January 2007. 
Near IR spectra were obtained of four jets from Class and I 
sources (see Section|2]i for one slit position at a given distance 
along the outflow, and orientated perpendicular to the direction 
of propagation. Observations were made (with position angle 
measured east of north) of the HH34 jet at 1" (400 AU) from 
the source (P.A.,/,,=77, 257°according to guide star availability 
on the observation date), and the H knot in the HH 1 1 1 jet located 
at 34" ( 14000 AU) from the source (P.A. sUt =T). Observations 
were also made of one knot in the approaching and receding jet 
from HH212, namely NK1 and SKI respectively, both located 
at 6" (2500 AU) from the source (P.A.. S ,„=1 16°). 

Using the long slit (49") and the long-blue camera configu- 
ration, in combination with the 110.5 1/mm grating and choos- 
ing a 2-pixel slit- width (0710), we achieve a spectral resolving 
power of R ~ 17 800. The corresponding velocity resolution is 
17kms _1 . Average seeing during the observations varied be- 
tween a/5 and a/8. 

We can use profile fitting to achieve an effective spectral res- 
olution which is beyond the actual resolution of the observations. 
Given a line profile which is intrinsically Gaussian, fitting allows 
an effective resolution which improves with climbing signal-to- 
noise. The one sigma error on the centroid is given by : 

resolution 

0~ centroid — , _ . „ \X) 

2.3548 snr 

where resolution is either spatial or spectral, de pending on the 
applic ation, and snr is the signal-to-noise ratio (IWhelan et all 
120081) . 

We chose filters G0504 and G0503 respectively, in order to 
detect the emission of [Fe II] 1 .64 fim and H2 2. 12/im. The stan- 
dard procedure for observations in the near IR was adopted. 
Eight integrations of 300 s were made, using the ABBA tech- 
nique, and then co-added. Two such exposures in each filter 
were obtained, in order to reach sufficient signal-to-noise in 
the jet borders. The data were calibrated and reduced accord- 
ing to standard procedures, using the tools within IRAF specif- 
ically designed for the reduction of GEMINI data. This in- 
cluded running tasks to correct for non-linearity (nsprepare) 
and spatial distortions (nsdist). Wavelength calibration was 
achieved using a dispersion scale derived from arc lines 
(nswavelength), which straightened and dispersion corrected 
the exposure (nstransform). All velocities are adjusted to 
the LSR reference frame (rvcorrect) and further adjustment 
was made for the LSR velocity of each system, Vlsr> TableQ] 
Additional overheads of standard star observations were not re- 
quired for a radial velocity analysis, and so data are not flux cal- 
ibrated. 



4. Data Analysis 

If protostellar jets are accelerated magneto-centrifugally then, 
with sufficient observing resolution and good signal-to-noise, it 
should be possible to detect the rotation signature of the flow. 
For good signal-to-noise, we can expect a typical accuracy in 
radial velocity of ~1.5kms _1 . In the case of a rotating jet, the 
gas on one side of the flow axis is expected to display a higher 
Doppler shift than on the other side. Therefore, measuring a gra- 
dient in the radial velocity profile in the direction perpendicu- 
lar to that of jet propagation can be interpreted as a signature 
of jet rotation. In such a case, the magnitude of the radial ve- 
locity differences (typically accurate to ~2kms -1 ) between the 
two sides of the flow axis allows derivation of a toroidal veloc- 
ity component. Observationally, this requires that the jet axis is 
identifiable, and that high spatial and spectral velocity resolu- 
tion is achieved to reveal the velocity profile as a function of 
distance from the jet axis. Clearly, this study is very demanding, 
and so obtaining useful observations from the ground is difficult. 
If the signal-to-noise is sufficiently high, however, we can reach 
beyond the nominal resolution via profile fitting of jet emission 
lines. 

As expected, the H2 emission was found to be spatially 
broad, and so is almost always spatially resolved (with the ex- 
ception of HH 34). Therefore, we could profile fit in the disper- 
sion direction to obtain a jet radial velocity profile transverse to 
the flow propagation, and thus determine the magnitude and di- 
rection of any implied jet toroidal velocity. First, the emission 
was binned and fitted with a Gaussian in the spatial direction, 
the peak of which was assumed to denote the jet axis. This was 
used as a reference point in determining differences in radial ve- 
locity between one side of the jet and the other. Each pixel row 
of the unbinned emission was then fitted with a Gaussian in the 
spectral direction, and the measured radial velocities at mirrored 
positions either side of this central position were subtracted. In 
a couple of cases where the emission profile was not entirely 
symmetric, cross-correlation was deemed more appropriate. In 
this way, we find how any differences in Doppler shift vary with 
distance from the jet axis. 

For the [Fell] emission, the spectral profile clearly traces 
both a high and low velocity component of jet material. Spatially, 
however, the emission is more confined to the jet axis and the 
seeing is never sufficiently clear for emission to be resolved im- 
plying the atomic jet width is at least <a/6. Thus, in our see- 
ing limited data, we cannot access the Doppler profile across the 
atomic component of the flow. 

5. Results 

5.1. HH34 

The HH 34 jet was clearly detected in both [Fe II] and H2 emis- 
sion. The observations are taken at only 1" (400 AU) above the 
disk-plane. Exposures for a given filter were not co-added, in 
order to preserve the information obtained in higher seeing con- 
ditions. FigureQ] shows the position-velocity plots for one ex- 
posure, observed in October. The data contained reflected con- 
tinuum emission with which the seeing at the time of obser- 
vation could be measured. Unfortunately, the jet width is un- 
resolved in both lines under the seeing of a/4 and a/6 re- 
spectively. The [Fell] higher velocity component (HVC) is 
travelling a t -96 ( + 2) km s , consisten t wit h measurements of 
iDavis et ail d2003l) : iTakami et alj (2006) and lGarcia Lopez et all 
(2008). Meanwhile, the Hi lower velocity component (LVC) is 
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Target 


RA 


Dec 




Vlt, 


' jet 


References 




(2000) 


(2000) 


(deg) 


(kms- 1 ) 


(deg) 




HH34 


05 35 29.8 


-06 26 58.0 


166 


8.7 


23-28 


1, 2, 3, 4 


HH111 H 


05 51 44.2 


+02 48 34.0 


97 


23.0 


10 


5,6 


HH212NK1 


05 43 51.547 


-01 02 48.0 


23 


1.7 


4 


7, 8 


HH212 SKI 


05 43 51.246 


-01 02 58.35 


23 


1.7 


4 


7, 8 



Table 1. Targets investigated in this paper, all located in Orion at 414pc dMenten et all 120071) . Our slit is placed perpendicular to 
the position angle (east of north) of the jet/outflow, PA/; olt ,. The systemic velocity is taken as the cloud LSR velo city, vi sr . The angle 
of inclina tion of the jet, ; ,>,, is give n with respect to the p lane of the sky . References - (1) Reipurth et al., 200 2t (DlAng lada et al 
fl995t (31lEisloffeletallll99l (4) iHeathcote etafl [l992l (51 iReipurfh L [l989l (6) iReipurfh et all Il992h (71 IClaussen et alT 1l998 
(8l lWiseman et alll200ll 
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Fig. 1. Position-velocity diagram of [Fe II] and H2 emission from the HH 34 approaching jet. Positive distances correspond to 
south-west. Both are spatially unresolved under the seeing of 074 and 076 respectively. 



travelling at -15 (+l)kms _1 . The H2 exposures show a decelera- 
tion from -15 to -10 kms 1 between October and December ob- 
servations. There is currently no available measure on the sense 
of disk rotation for this system. 

5.2. HH 111 - H 

While the HH 1 1 1 jet is clearly detected in both [Fe II] and H 2 
emission, Figure|2] the signal-to-noise in H2 was extremely poor 
in spite of the long integration. There was no reflected contin- 
uum emission in either line from which to determine the see- 
ing during observations, and so we rely on the weather report 
at the time. The [Fell] HVC is travelling at -107(±l)kms _1 , 
while the H2 LVC is travelling at -39 (+2) kms 1 , in line with 
iDavis et a il (I2001bh . Disk observations show that the northern 
side is blue-shifted w ith respect to the southern side ( Yan g et all 
ll997tlLee et alll2009l) . We therefore could expect to detect pos- 
itive values of radial velocity differences, Av ra ^. Our results are 
all within error bars about zero, and so we do not detect a clear 
Doppler gradient. A mere hint of a gradient may be apparent, 
but these are in the opposite direction (i.e. negative values of dif- 
ferences in radial velocity, Av ra d) to that expected based on the 
sense of disk rotation measurements. Note that we are examining 
this jet knot at 33" from the source and so environmental influ- 
ences are more likely come into play in significantly disrupting 
the flow kinematics. Finally, in the position- velocity diagram, we 
see that a second emission peak in H2 was observed (offset spa- 
tially by 0"2). The emission is much fainter and appears to trace 
an offset HVC in H2, with Wlsr measured as -110(±3)kms _1 . 
A similar faint HVC is some t imes s een in echelle spectra of bow 
shocks, see e.g. lSmifh et alj d2003l) (on HH 7). In this case, the 
HVC is attributed to a fast-moving mach disk. 



5.3. HH212 - NK1 

There is a clear detection in both [Fell] and H2 emission in 
the NK1 knot from the HH212 outflow of its deeply embed- 
ded Class source, Figure[3] [Fe II] emission is notably fainter 
with about a third of the signal-to-noise of H2. The [Fe II] HVC 
is travelling at -27 (±l)kms _1 , while the H2 LVC is tr avel- 
ling at -8(+l)kms -1 , consistant with IDavis et al.1 (l2000l) and 
ICorreia et alj (12009). The jet width is unresolved in [Fell] emis- 
sion, but resolved in H2, and shows a slight hint of a Doppler 
gradient. Although nearly all Av rai / data points are within error 
bars about zero, there is a systematic trend such that the dif- 
ferences in radial velocity are consistently negative rather than 
randomly scattered about zero. The Doppler analysis suggests 
that the south-east side of the flow to be redder than the north- 
west side. However, for NK1, this spatial intensity distribution 
is asymmetric (as seen in t he position-velocity plot), and as 
also seen in the H2 study by ICorreia et all d2009l) . Recall from 
Section|4]that the emission was binned and fitted with a Gaussian 
in the spatial direction to determine the jet axis, and that this was 
used as a reference point in determining differences in radial ve- 
locity between one side of the jet and the other. Therefore, the 
spatial asymmetry introduces inaccuracies in determining the lo- 
cation of the jet axis, and hence the velocity differences across 
the jet. 

5.4. HH212 - SK1 

The receding jet knot SKI from HH212 was clearly detected in 
both [Fe II] and H2, Figure|4] Again, the jet width in [Fe II] emis- 
sion is unresol ved. The H2 LVC is travelling at +6(+l)kms~', 
consistant with IDavis et all (I2000I) . For H2, the values of Av ra d 
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Fig. 2. Top: position-velocity plots for the H knot of the HH 1 1 1 jet in [Fe II] and H2 emission, where positive distances correspond 
roughly to north; Bottom: radial velocity profile across the jet in the spatially resolved H2 emission, and a plot of the differences in 
the Doppler shift, Av ra j, between one side of the jet knot and the other. 




Fig. 3. Top: position-velocity plots for the approaching lobe of the Class bipolar jet HH 212, NK1 knot, in [Fe II] and H2 emission, 
where positive distances correspond to south-east; Bottom: radial velocity profile across the jet in the spatially resolved H 2 emission, 
and a plot of the differences in the Doppler shift, Av ra ^, between one side of the jet knot and the other. However, caution must be 
taken in its interpretation given that the emission intensity is asymmetric in the position direction. 
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are consistently negative, although again we see a spatial asym- 
metry in the emission (see position-velocity plot). The direction 
of the Doppler shift reveals the south-east side of the flow to 
be redder than the northwest side. This matches the hint of a 
gradient seen in our NK1 dataset (Section l531 ). It also ma tches 
the direction of the H2 measurements of iDavis et al. (2000) and 
ICorreia et al.l (120091) . and the direction of th e Doppler gradient 
reported for the disk dWiseman et all 12001 1) . However, caution 
must be taken in the interpretation of the H2 gradient for both 
NK1 and SKI. The emission intensity is asymmetric in the posi- 
tion direction in both lobes, bu t on opposite sides o f the jet axis, 
as also seen in the H2 study bv lCorreia et al.l d2009l) . 

Overall, the results for HH 212 are in agreement with those 
obtained bv lDavis et al.l (2000) for SKI with differen t teles copes 
and at different epochs (7 years apart). IDavis et all (120001) and 
iLee et al] (120081) both report a high positive Doppler slope in 
SKI with respect to NK1, and we find the same trend. Therefore, 
these signatures are real and persist over time. 

Lastly, note that the [Fell] HVC is measured as having a 
radial velocity of -25 (±l)kms _1 . This is highly unusual given 
that this is the receding jet lobe. To ensure the measurement is 
correct, we have checked the offsets of the slits in the raw data 
headers; the LSR velocity correction with respect to the obser- 
vation dates; the position of the emission line with respect to the 
OH sky lines in the original raw data images; the wavelength 
calibration procedure; and the velocity of the H2 line with re- 
gard to the literature. No errors were found in the acquisition or 
calibration of the data, and all other measurements are consistent 
with previously reported values in the literature. Therefore, this 
velocity measurement appears to be real but we cannot find an 
explanation for it based on these data. 

5.5. Sources of Errors 

Difficulty in determining transverse Doppler profiles arises if the 
slit is unevenly illuminated, e.g. if the target is not well-centered 
in the slit. Using a slit orientated perpendicular to the jet prop- 
agation direction (as in our observations) ensures that uneven 
slit illumination is a second order effect, and hence the prob- 
lem o f the so-called slit-effect is negligible dChrvsostomou et all 
120081) . This was not the case in previous studies using a slit par- 
allel to the flow (i.e. iBacciotti etaTl 120021; IWoitas et ail 120051) . 
and so the subtraction of the uneven illumination effect added 
an element of complexity to the data reduction which we have 
avoided here. 



6. Discussion 

We have detected gradients in the transverse Doppler-shift pro- 
file in resolved H2 jet targets. We now examine whether or not 
we are justified in interpreting these gradients as a rotation of 
the flow, and thus whether we can claim to support the magneto- 
centrifugal mechanism of ejection. We first investigate whether 
the implied toroidal velocity and angular momentum flux are re- 
alistic, in the context of jet and disk rotation measurements re- 
ported in the literature. We then examine the statistical founda- 
tions for such a claim. 

The implied jet toroidal velocity may be calculated from the 
measured radial velocity differences. The implied toroidal veloc- 
ity is = (Av rac i/2)/ cos ij et , where ij et is the inclination angle 
with respect to the plane of the sky, and Av rm / is the radial veloc- 
ity difference between two sides of the jet. We see a systematic 
radial velocity difference of ~3kms _I for the HH 212 jet. The 



implied toroidal v elocity is then ~ 1 .5 km s . T his is in line with 
values derived by Chrvsosto mou et a"D (120081) of 1 to 7 kms 1 
measured for two Class I sources examined in H2 emission. 
These authors demonstrated the velocity signature to be consis- 
tent with a simple jet rotation model. The value is also in the 
same range as t hose m easured in H2 emission for HH 212 SKI 
by IDavis et all (l2000l) . in which the authors point out that the 
expected jet rotation speed bas ed on the disk rotation measure- 
ments of lWiseman et al.l(l2001l) is ~2kms~', thus demonstrating 
that our measurements are consistent with a jet rotation interpre- 
tation. It is also clear from the modeling of Corr eia et al. 1 (120091) 
that the toroidal velocity may be confused with other kinematic 
signatures, and so our calculations are broad indications. 

Comparing with more evolved Class II TTauri jets, our 
toroidal velocity is substantially lower than those measured 
in optical and near-UV emission close to the launch point , 
which yielded typically 10 to 20 kms~ 1 dBacciotti et all (2002; 
IWoitas et all 120051: ICoffev et all 12004 2007). These results 
originate from the resolved atomic emission which is more colli - 
mated, and so would be expected to yield higher toroidal velocity 
values with respect to the molecular component. Unfortunately, 
without resolving the atomic component in our Class 0/1 sample, 
it is not possible to make a direct comparison between the two 
classes, in order to understand the evolution of angular momen- 
tum extraction over the age of the source. 

Nevertheless, we attempt to gain a rough estimate. 
The angular momentum extraction may be appro ximated as 
Ljet ~ rv^Mjef. For the molecular knot HH212 SKl. lDavis et al.l 
(120001) report M jet ~ 4.2xl0~ 8 Moyr 1 . This is in the 
range of 2 to 5xl0~ 8 M H yr 1 reported for Class I outflows 
(lAntoniucci et all 120 08). We adopt V0=1.5kms _1 at an average 
distance from the jet axis of Of! 4. Although the jet is resolved, we 
do not see reflected continuum emission which we could use in a 
PSF deconvolution to determine the true jet radius. These values 
give Ljet ~ lxl0~ 5 M yr _1 AUkms -1 . This is in line with that 
for HH26 of 2xl0~ 5 Moyr 1 AUkms 1 derived at 1-2" from 
the disk plane, for a jet radius of 0744. However, our result im- 
plies a specific angular momentum whi ch is 3 times higher than 
that of the SiO knot SS lLee et aD (f2008). which lies closer to the 
star at 2". The difference mainly arises from a smaller adopted 
SiO jet radius of 071. 

Roughly comparing Class 0/1 values with those for Class II 
jets, we take the example of TTauri systems CWTau, RW Aur 
and DG Tau, which have been studied c lose to the je t foot point, 
i.e. a/5 (70 AU) from the disk-plane. ICoffev et all (120081) re- 
port angular momentum fluxes in one lobe from each system 
as 1, 3 and 13 xl0~ 6 M yr _1 AUkms -1 respectively for the 
atomic component. We find up to an order of magnitude dif- 
ference between the angular momentum extraction of Class 0/1 
versus Class II sources. This is a rough indication of the mag- 
nitude of the decrease in angular momentum extraction as the 
young stars evolve. Indeed, the decrease in mass accretion flux 
suppo rts this: a co mparison between Clas s II dGullbring et all 
1998) and Class I ( Antoniuc ci et all [2008) shows the mass ac- 
cretion rates for the Class I are between one and two orders of 
magnitudes larger than those of Class II. 

It seems possible that we are observing a jet rotation signa- 
ture in these data. We must also consider whether it is probable. 
Does the gradient in the radial velocity profile in the transverse 
direction represent a rotation of the jet, or could it be a signature 
of an asymmetric bowshock, precession of the flow, or indeed a 
combination of effects. 
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Fig. 4. Top: position-velocity plots for the approaching lobe of the Class bipolar jet HH 212, SKI knot, in [Fe II] and H2 emission, 
where positive distances correspond to south-east; Bottom: radial velocity profile across the jet in the spatially resolved H2 emission, 
and a plot of the differences in the Doppler shift, Av ra d, between one side of the jet knot and the other. However, caution must be 
taken in its interpretation given that the emission intensity is asymmetric in the position direction. 



Steady-state MHD jet models imply that the rotation of the 
flow should persist to large distances, as a necessary outcome 
of magnetic flux conservation. Nevertheless, an important ob- 
servational precaution is to examine the flow close to the launch 
point. In this way, we try to ensure that environmental factors 
do not significantly disrupt the inherent jet kinematics. This is 
not always observationally possible due to high opacity close 
to the star, especially when dealing with younger more embed- 
ded flows as is the case here. Hence the need to observe these 
younger outflows at the position of the brightest knots, which 
are not necessarily close to the launch region. This observational 
obstacle we must examine the overall context of our results. 

Support for a rotation argument lies in important consistency 
checks. For example, if the jet is rotating as a result of the launch 
mechanism, we should consistently measure gradients in the jet 
Doppler profile transverse to the flow direction for many targets. 
Also, the direction of the implied jet toroidal velocity should 
match measurements for the direction of the disk rotation within 
the same system, since the jet is supposedly extracting its angu- 
lar momentum. Likewise, the implied sense of rotation should be 
the same in both lobes of a bipolar jet and should persist along 
the flow. 

From our survey, it appears that transverse Doppler gradients 
may be consistently measured in Class 0/1 jets in the near IR H2 
e mission, as we set out to esta blish following the initial findings 
of IChr vsostom ou et al.l d2008l) . Furthermore , we provide confir - 
mation of a previous detection for HH 212 Ibavis et al (2000), 
thus also demonstrating a persistence of the transverse Doppler 
gradient over time (from 1999 to 2006). Together with the fact 
that gradients are also consistently measured in Class II jets in 
the optical and near UV, we have an encouraging case so far. 



Next we consider Doppler gradients of jet and disk. 
Agreement is compulsory to validate a jet rotation interpreta- 
tion. For HH 2 12, we confirm an agreement in sense, as also re - 
ported by both lDavis etal] d2000t) for SKI and lLee et al1d2 008). 
Howev er, they do not agree with the NK1 results of Davis et akl 
(120001) or the SiO results of ICodella etal] d2007l) . Furthermore, 
given the spatial asymmetries we find in this case, our results 
must be treated with caution. We cannot determine a clear 
gradi ent in the jet for HH 111 in order to compare with the 
disk ( iYang et all Il997t iLee et all 120091) . but the data suggest 
a gradient which opposes that of the di sk. Other cases exam- 
ined are the TTa uri systems DG Tau dB acciotti et al., 2002; 
i Testi et all l2002h. RWAur (ICoffey et all 12004 ICabrit et al 
2006), CWTau icoffeyet all 120071) and HH30 (ICoffev et al 
l2007t iPetv et all l2006). Agreement has also been confirmed 
for DGTau, in which jet observations were conducted close 
to the launch point and where the transverse gradient persists 
for 100 AU. Disconcertingly, RWAur shows a disagreement in 
sense. This may be a complex triple system and hence must 
be further investigated to confirm that no other influences are 
coming into play. For CWTau, agreement has been found in 
preliminary disk results (C. Dougados, private communication). 
For HH30, the jet rotation sense was deemed inconclusive. 
It has since been revealed that the HH 30 jet is in fact wig- 
gling ( Anglad a et all 20071). due to th e orbital motion of the 
binary source (Guillote au et all 120081) . Lastly , agreement has 
been recently found in both HH 21 1 and CB26 ( iLee et all 120071; 
lLaunhardt et al.Ll2.Q09t) . Overall, of the 8 systems for which both 
jet and disk gradients are studied at this early stage of our work, 
four show clear agreement and 1 shows clear disagreement. The 
statistics are as yet too low to be significant. 
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Agreement between Doppler gradients in both lobes of a 
bipolar jet is also compulsory to validate a jet rotation inter- 
pretation. We measure an agreement in the gradient between 
the two lobes of the HH212 bipolar flow. However, we again 
sound a note of caution based on the spatial asymmetry and 
the fact that we observe far from the source. This supports the 
agreement also f ound in SiO measurements in the two lobes 
dLee et al.L 1200 8), although the same agreement was not found 



Davis et alj d2000). Only two other bipolar flows have been 



examined. They are from the T Tauri systems Th 28 and RW Aur 
dCoffev et all £004), and in both cases agreement was found. 
Furthermore, in two cases were these bipolar jets have been ob- 
served with the slit parallel to the flow, the gradient has persisted 
along the jet over a distance of 90 AU in the same direction in 
both lobes. Such a scenario is not likely to arise from a signature 
of asymmetric shocking, and certainly not from jet precession 
where the gradient should be opposite in the two lobes. In other 
words, in the case of precession, for several slits parallel to the 
jet (mimiking an IFU), the poloidal velocity peak is offset with 
respect to the central slit in opposite directions for each lobe, 
and hence the radial velocity profile shows a gradient opposite 
in direction in each lobe. In this case, also, there is a periodic 
change in the direction of the gradient as a function of distance 
along each jet lobe. Such changes are not found in these parallel 
slit studies. 

Although statistics are limited, we are building our way to- 
wards a statistical argument to support the fact that we are in- 
deed observing jet rotation. Such information is critical in fi- 
nally providing observational confirmation for the widely ac- 
cepted but untested centrifugal MHD wind launching mecha- 
nism. Obviously, any evidence which suggests protostellar jets 
are launched centrifugally would, of course, in turn support the 
idea that the same mechanism is at work in their larger-scale 
brethren, i.e. the AGN jets. Indeed, the latest observations in- 
dicate jet rotation is also detectable in active galactic nuclei 
dYoung et al.L 120071) . Together with observations indicating jet 
rotation in protostellar jets, this provides support for the univer- 
sality of the theory of magneto-centrifugal ejection. 

7. Conclusions 

We have conducted a ground-based survey in near IR lines of 
four jets from Class 0/1 sources to search for signatures of jet ro- 
tation. These embedded sources make the use of adaptive optics 
very difficult, due to the absence of a nearby optical guide star. 
Thus, our results are derived from seeing-limited observations 
and rely on profile fitting in the spectral direction. For HH 34, 
we find emission in [Fe II] and H2 is detected but spatially un- 
resolved with seeing of 0'.'4-0'.'6 across the jet within 1" of the 
driving source. For HH 1 1 1-H, we find we cannot resolve a clear 
gradient across the jet, but perhaps glimpse only a hint of one, 
which is opposite to the direction of the disk. Similarily, in the 
NK1 knot of the HH 212 bipolar flow, we see a hint of a gradient 
but this time in the same direction as the disk. Lastly, we detect a 
gradient of typically 2-5 kms 1 in SKI knot of HH212. We find 
agreement in the sense of the gradient in both HH212 lobes, 
as would be expected if the measurements are of j et rotation. 
The re sult confirms an earlier detection in SKI bv iDavis et alJ 
d2000l) an dlCorreia et alJ (|2009l). and also matches the disk rota- 
tion sense I Wiseman et all d2001l) . Furthermore, it demonstrates 
the persistence of the gradient over a time frame of seven years. 

For HH212, we can estimate a possible implied toroidal 
velocities of 1.5kms _1 , previously sh own to be consis- 
tent with a simple jet rotation model dChrvsostomou et all 



120081) . and angular momentum flux estimates of Lj et ~ 1x10 
M Q yr _1 AUkms -1 . A direct comparison cannot be made be- 
tween Class 0/1 and Class II angular momentum extraction, since 
[Fell] is unresolved in our data. We are therefore prevented 
from determining any particular trend over evolutionary time. 
However, for the HH212 and HH111 knots, which were ob- 
served very far from the driving source, we must consider that 
the gradients are less likely to be unaffected by external factors 
such as asymmetric shocking. 

Our analysis illustrates the observational difficulties in con- 
ducting this study, as it pushes instrumentation to the limit in its 
demand for a combination of high spatial and spectral resolu- 
tion, as well as good sensitivity. Overall, it is clear that in order 
to safely interpret Doppler gradients as signatures of jet rotation, 
there is a need for improved statistics via high resolution multi- 
wavelength jet observations close to the source, and comparison 
with the associated disk rotation sense. 
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